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1 Introduction 

Recent studies[|ll-|0| reveal an unexpected connection between a certain category of or- 
dinary differential equations (ODE) and integrable models (IM) with quantum group 
symmetry. We call this ODE/IM correspondence Q The success inherits the fruitful re- 
sults from the exact WKB analysis [3-|[T6[| and progress in the study of integrable structure 



m-m- 

The aim of the present talk is two-fold, the survey on the 2nd order ODE case, and 
the brief preview on the generalization to higher order cases. 

Firstly, we will review the results on the 2nd order differential equations, the 
Schrodinger equation with a polynomial potential, 

+ X y = hy. 



dx^ 



We will mainly follow the argument in but employ some simplifications and add some 
materials. Motivated by the success of the exact WKB method, we regard the coordinate 
X as a complex variable. The complex x plane is conveniently divided into sectors. See 
section ^j. Each sector possesses two (= the order of the equation ) linear independent 
solutions. We call them a fundamental set of solutions (FSS). The relations among FSS 
of different sectors are of our interest. To be precise, we would like to evaluate the Stokes 
multiplier which characterizes the connection rule. In this view point, it is natural to regard 
that the problem consists of two coupled equations, the original differential equation and 
the difference equation for the Stokes multiplier. It will be then shown that fundamental 
objects in the study of the solvable models, e.g., Baxter's Q— operator, fusion hierarchy 
of transfer matrices based on Ug{A\ ) and their functional relations naturally come into 
play in the analyses on ODE. Especially, (unfused) transfer matrix is identified with the 
Stokes multiplier. Reflecting the ODE/IM correspondence, the Stokes multiplier has two 
representations, the Wronskian representation, which arises from ODE, and the DVF 
representation, originated from IM. They both play a role in generalizing the results in 
the second part of the talk. 

Physically, the Stokes multiplier may be less interesting. Rather, the quantity of im- 
portance is the spectral determinant, D{E) = det(?^ — E) or eigenvalues themselves. 
Remarkably, D{E) also belongs to the fusion hierarchy. Thus the result provides a unified 
view of Stokes multipliers and spectral determinant. 

In the second part, a generalization to higher order ODE will be addressed, 

dJ'+^y . 

+ X y = Ey. 



dx^' 

In Q, functional relations were derived among Stokes multipliers and their general- 
izations. These are identical to functional relations among transfer matrices of solvable 
models with Ug{An ) symmetry, which generalizes the observation for n = 1. The rela- 
tions were evolved by use of the machinery in the solvable models, a quantum analogue 

^Following R. Tateo. 



of the Jacobi-Trudi formula. Here we will give an alternative, much simpler derivation, 
resulting from the Wronskian representation of the Stokes multipliers. We also note that 
the possible connection of the relations to discrete soliton equations (the Hirota-Miwa 
equation) [pl , ^ . 

The parallelism to IM will be further exploited. The eigenvalues of transfer matrices 
possess a universal structure called the dressed vacuum form (DVF). The universality 
has a deep origin in analyticity of their expressions under Bethe ansatz equations and 
the Yang-Baxter integrability. We will show that Stokes multipliers also assume the same 
DVF. 

The paper is organized as follows. The asymptotic form for n + 1— th ODE will be 
discussed in section |[ Several notations and symbols, such as sectors, Stokes matrices, 
are introduced for n general. In the next three sections, we restrict ourselves to the n = 
1 case. The recursion relations and functional relations for the Stokes multiplier and 
its generalizations are derived in section ^. Under certain assumptions, one transforms 
the algebraic relations to a set of integral equations modulo one unknown parameter. 
Remarkably, the integral equations take identical forms to thermodynamic Bethe ansatz 
equations. We shall discuss the DVF representation for the Stokes multiplier in section 
^. The spectral problem is addressed in section ^. Utilizing the previous results, spectral 
determinants are identified and the one missing parameter in section ^ will be determined. 
The extension to arbitrary n is the topics of sections ^ and |^. We conclude the paper with 
a brief summary and discussion in section B 



While preparing the manuscript, I find the preprint p3| appearing on e-print. The 
content of the paper largely overlaps with the second part of the present manuscript. They 
actually treated a more general set of ODE but without the argument of the functional 
relations. 



2 Asymptotic Expansion, FFS and Stokes multipli- 
ers 



The details of the present section can be found in [|J, |2^, |2^ . 



We first discuss the asymptotic behavior of a slightly generalized differential equation, 

<9"+^2/ + {-lTP{x)y = (1) 






i=o 

(2) 

where Qj are complex numbers and ai = 1. Note that the factor (—1)" is not essential. It 
can be adsorbed into re-definition of the angle of x. For later convenience, we will include 
this factor throughout out this report. 

Now that X = cxD is an irregular singular point of the equation, analytic properties 
of the solutions are different for different angle regions in complex x plane. Let Sk be a 



region in the plane satisfying 

77 



I argx + k9\ < 



i + n + 1 



for X G 5fc, where 6 = j^ri- We first analyze the asymptotic behavior of a subdominant 
solution in Sq. Following [^ ^], we define hh{h = 1, 2, ■ ■ ■ ) by the relation, 



(1 + J] afcX-'=)i/("+i) = 1 + 1] ^hX-\ 
k=i h=i 

A key function E{x, a) is defined by bh, 

/ he 
h=i 



^ + 1 ,(^+n+l)/(n+l) _^ V^ '^^ ^£/(n+l)+l-fe 



_^,. r. .,.;/,.'.-;+ X j.^ 



e + n + 1 ^_£__/, + i 

where hi = N ioi i = N{n + 1) — j, (j = 1, ■ ■ ■ , n). a stands for (ai, 02, ■ ■ ■ , a^) 
In addition, we introduce an exponent ue by 

f, for £^Omodn + l ,. 

^^~ f + (n + l)6?,,+i, for £ = mod n+1 . ^^ 

Theorem 1. In So, there exists a subdominant solution to ^ y{x,a.) which has the 
asymptotic behavior, 



y{x, a) ~ c-'^x~'^i/in+i)^-E{x,^)_ ^4) 

A normalization factor C is introduced for convenience in the later discussion, 

(7"+^ := exp( i) TT (w^-w'), w := exp( i). 

2 -'--'- n+1 

0<i<j<n 

As argued in [^], the range of the validity of the asymptotic form is wider if one forgets 
the subdominance. Explicitly, it is valid for |argx| < i^^+i '^- 

The intriguing feature in the differential equation (|T]) is a certain symmetry in rotating 
X plane. 

Theorem 2. if y{x,a) is the prescribed solution, then 

2/fc(x,a):=l/(xg-^GW(a))g"'=/2 

is also a solution to ^. 



The parameter q signifies exp(z6') = exp(z ^_^^^_^-^ ). The operation G^^\a) is defined by 

G(fc)(a) = G(G('=-i)(a)),fc>2andG(a) = (ai/g, a2/g^ ■ ■ ■ a^/gO- 

From now on, we restrict our discussion to a single potential term case , 

P(x) = x^ + a£, a^ = A"+^ 

One immediately verifies that &/i^+i = and thus v^ = ni/2 for i > n + 1. Under the 
operation of G, G{ai) = aeq~^ = a^g""^^. In term of A, the action of G is simply given by 
G'('=)(A) = Ag''. Consequently, yk = q'^^/'^y{xq~^,\q'^). 

A set of fundamental solutions (FSS) in Sk is formed by by {yk,yk+i, ■ ■ ■ ,yk+n)- We 
introduce a. (n + 1) x (n + 1) matrix ^kix) 



$fc(x) := 



/ Vk, Vk+i, 
dyk, dyk+i, 

yd^'yk, d'^yk+i 



Vk+n \ 

dyk+n 



dyk+nj 



(5) 



We denote the Wronskian, the determinant of ^k{x), by Wk- Note that the above 
asymptotic expansion is valid for yk+j, (j = 0, ■ ■ ■ , n) in the common sector Sk+i/2^<Sk-i/2- 
As Wk is constant in x, one easily checks the linear independence of these solutions by 
using the asymptotic expansion (^ at the sector. Due to the present normalization oi yk, 
we have Wk = 1. 

A Stokes matrix 5*^ connects FFS of Sk and Sk+i 



$fe+l(,T) = ^k{x)Sk. 
The linear independence of solutions demands 5*^ in the following form. 



(6) 



/ rfHAg'^), 1, 0, 0, ••-, OA 

rfHAg'^), 0, 1, 0, ■■■, 

Sk = ': : 

r^^^lAg'^), 0, 0, 0, ■■■, 1 

W^'\M'), 0, 0, 0, ■■•, 0/ 

We call elements r Stokes multipliers. 

By the Cramer's formula, one represents t^^ (Xq'^) as 



(7) 



-iJ)f\^k 



Tl"{Xq'') = det 



Vk+l, yk+2, 



5"|/fc+i, 9"|/fc+2, 



Vk, 

d'^yk, 




that is, {i/k, dt/k, • • • , d^-yk) is inserted in the j— th column in the denominator. Evidently 
ri^+'\Xq'^) = (-l)nw,_,,/W, = i-ir. 

The above representation (|]) of Stokes multipliers will be referred to as the Wronskian 
representation. 

Determinants of such structure will be hereafter abbreviated to, by specifying only the 
first row, [yk+i, yk+2, ■■■ ,yk,--- , Vk+n+i]- Generally, 

Vil ) 1/J2' ■ ■ ■ ; Vin 

[yiii yi2i ' ' ' 1 Uini -^ : : • 

\a"-iy,„ d^-'y,,, ■■■, d^'''y.,J 

We have prepared materials needed for study on general n. In the next few sections, 
however, we confine ourselves to the n = 1 case. There are two reasons for the separated 
argument. First, only for n = 1 case, we have a clear bridge between the connection 
problem and the spectral problem. Second, the second order ODE may be the most 
relevant to physics. 

3 Fusion Stokes matrices for the 2nd order ODE 

The ingenious idea in lies in the introduction of the generalized (or fusion) Stokes 
matrices connecting the second neighboring sectors, the third neighboring sectors, and so 
on. We denote by Sj^ the fusion Stokes matrices connecting two FSS, $fc and ^k+j 

$, = <^k+,Si^\ 
Obviously, the recursion relation holds, 

Sl^^ = Sii\'^S^^\ (9) 

Theorem 3. Sf^ has an expression 

where we adopt Tq (A) = l,rl/(A) = 0. Thanks to the condition ye+2+k = —Vk, t) (A) = 

— rj^2('^) = 1 ^'^^ 'i£+i('^) ~ '-'■ Naturally, r- , (j > 2) are referred to as the generalized 
Stokes multipliers. Due to ([^ they satisfy relations, 

r«(gA)r«(A)=r«,(A) + rji\(g2A). (10) 



example 



For £ = 1, 

~--y + xy = X y, 
ax 

it is well known that the eigenfunction is given by Airy function y = Ai(a;). The above 

connection rule then fixes the Stokes multiplier for Airy function rj = 1. 

Let T^^\X) = TJ^^(Ag(^+i)/2)^ Que can then prove 

^(1)^^^1/2)^(1) (^^-1/2) ^ 1 ^ i;« (A)T,(1) (A). (11) 

using ([T0|) and the mathematical induction. These functional relations exactly coincide 
with those among fusion transfer matrices ofUq{A[ ). In the latter context, the suffix j 
specifies the spin j/2 assigned to the auxiliary space. They are the closed set of equations 
among finitely many unknown functions T , (j = 0,1,- ■ ■ ,i). Thus they may be of 
significance in the estimation of the quantity of our original interest, t{ (A). Actually, 
with additional assumptions on the analyticity and asymptotic behavior of r- (A), one 

can fix r} (A) via coupled nonlinear integral equations resulting from ( |ll]). To see this, we 
conveniently put Yj(\) = Tj_|_i(A)T,_i(A) and A = e^/'^^^^-'^. Now the functional relations 
read 

y,(^ + ^^)F,(^-zp = (i + Vi(^))(i + r,+i(0)), J = 1,2, ■■-,£- 1 (12) 

where h = i. Note that Yo = Ye = 0. 
Assumption 

(logy,(^))' ((log(l + Yj(6'))') are analytic, nonzero and have constant asymptotic behavior 
(ANZC) in the strips Im^ E [-f , f], Im^ E [-0+,0+] respectively. 
The validity of this assumption will be discussed in section ^ 

Once this is granted, one immediately derives from ([12|) p, ^ 0, 



^ /i-i 



e .(61) = rrijr expO - --} (j)j,k * Lk{6) (13) 



k=l 



where ruj = sm{7ij/h)/ sm{n/h), Yj{9) = exp{ej{9)) and Lj{9) = log(l + l/Yj{9)). 
The asterisk denotes the convolution, A * B{9) = J A{9 - 9')B{9')d9'. 

This type of coupled integral equations is known as thermodynamic Bethe ansatz 
equation(TBA). One finds them in various branches of IM, e.g., the thermodynamics of 
ID spin chains or the perturbation theory of CFT. The kernel, (pj^k, is related to the two 
particle S-matrix Sj^k of quantum field theory based on Ah-i by (t)a,h{9) = —ide log Sj^k{9) 



and 



min(jf,fc)— 2 

Sj,k{o) = n {ij-^i+2j+i}, 

3=0 
r T r .^r .s / n Smh{e/2 + i7l/2h) 

Theorem 4. The set of equations j^lBj ) fixes Yj{6) for a given r. 

To determine the factor r, it needs an independent ingredient from the spectral theory. 
We will come back to this point in section ^. 

We have a remark. For the later use in the spectral problem, we have introduced gen- 
eralized Stokes matrices and derived functional relations (|Tl|) from the obvious recursion 
relation (^. They can be also easily extracted from the following Wronskian representa- 
tion of rj^^(E), 



r^}\E) = det (yj^ yj-^\. (14) 

For n > 1, the situation is different. The generalized Stokes matrices can be defined 
similarly. Their elements, however, do not contain nice generalization of r|" s. The formal 



definition of the Wronskian type like (|Tj) still works efficiently. See the discussion in 
section ^. 

4 Dressed Vacuum Forms of Stokes multipliers 

As shown in the previous section, the Stokes multipliers share same functional relations 
with the transfer matrices of IM. Below we will discuss if this correspondence carries 
forward. 

The eigenvalues of the transfer matrices in solvable models exhibit a universal structure 
often referred to as the dressed vacuum form (DVF). We shall explain DVF for the simplest 
the A\ case with the dimension of the auxiliary space being 2. 

Obviously, the highest weight state (= vacuum) is the trivial eigenstate of the transfer 
matrix. Its eigenvalue consists of two terms, reflecting the dimensionality of the auxiliary 
space. Each of them is given by the simple product of the local weights which is termed 
as the vacuum expectation value, 

Tvacuum(A)=/l(A) + /2(A). 

This expression must be modified for general eigenvalues. The quantum inverse scattering 
method yields the exact expression. The result tells that Tvacuum(A) must be modified 



by "dressing " the vacuum expectation values with ratios of Baxter's Q operator (or its 
eigenvalue) which commutes with T, [T, Q] = 0, 



,-1^ 



,(..(„., ,„|^,,,„«^. (15) 

The fact that the eigenvalue must be pole free results the famous Bethe ansatz equation 
(BAE), 

MX,) _ Q{X,q-') 
U\,) Q{\q) 

where Q{Xj) = 0. This kind of representation is called DVF. 

Clearly, eq( |15D has an interpretation as the second order difference equation (Baxter's 
T-Q relation), 

Ti'\X)Q{X) = MX)Q{Xq) + f2{X)Q{Xq-'). 

Thus two independent solutions exist which we call Q±. 



In p7| , pH] , it is shown that DVF is universal for models based on general Uq^g) under 
certain assumptions. The key ingredient in the argument is the analyticity under BAE. 
Thus we may conclude that DVF embodies the BAE or Yang-Baxter integrable structure. 

Now we turn to the n = 1, /c = of (||). The Stokes multiplier r} ^(A) is given by y's 
in two manners, 



(1)..^ 1/0 + 1/2 , _ I/O + 1/2 I ..^^ 

2/1 Vi 

Originally, the rhs can be evaluated at any x yielding the same r| (A). We adopt a 
convention to enumerate them at the origin for the later convenience. See section ^ 

By comparison of (|15[) and ([T6|) and the identification, r} (A) = T^ (Ag), made after 
(p!0|) , we deduce yj oc Q^{Xq^) and y'j oc Q+{Xq^). Precisely, the argument in the next 
section concludes yj = q^^'^Q-{Xq^), and y'j = q~^^'^Q+{Xq^). 

The linear independence of FSS implies that r{ (A) is pole-free. On the other hand, 
y{0, A) can generally be zero for some A = A^. Thus we have BAE for Stokes multipliers. 

It is interesting that dy/dx, which is by no means a solution to the original ODE, 
now appears as the second " solution" to the difference equation. This issue will be further 
pursued in a later section. 

The coincidence is not only for the spin 1/2 case, but also for cases of arbitrary 
spins. This can be easily seen as they share the same initial condition and the functional 
relations. One can also verify this directly using the Wronskian representation. For this 
we rewrite the condition Wk = 1 in the form, 

/ I -, ^ Vk+l Vk 1 



2/fe+i yk ykyk+i 



With use of this, one obtains 



yj+i yo ^ Vj+i-k Vj-k 

^ 1 



= yoyj+i2_^ , 

^ yj^kVj-k+i 

which coincides with the known expression for the transfer matrix. Actually the dis- 
cussion like above has been firstly found as the operator identity under the name of the 
quantum Wronskian form. We follow the discussion in [0, ^ for reproducing the DVF 
for the spin j/2 case. 

Before closing the section, we present simplest examples (£ = 1,2 ) where explicit 
solutions are available by elementary functions [|I], ^ |^, |T^, |I3[. We shall use E instead of 



X {E = A^) and adopt same symbols, y, r etc as the function of E. 



The case 



This is a well known example in quantum mechanics. The wave function y is given by the 
Airy function. 

On the other hand, for £ = 1, we have r} = 1 (theorem 3). Thus T — Q relation 
simplifies. 



Q^E) = q'^'^Q_{Eq-^) + q^/^Q_{Eq^) 

Q+{E) = q'''Q^{Eq~') + q'^l'Q^{Eq'). (17) 

These relations coincide with the 3-solution dependence relation for the Airy 



function ||12||, 



q-^M{q-^E)+Ki{E) + qki{qE) = 
q-^Bi{q-^E) + Bi{E) + q'^Bi{qE) = 

where Bi(a;) := ^'^^ . 

To check this, we use g^ = 1 in the arguments in the first of ([T7|), 

Q4E) = q-^'^Q-{Eq) + q^/^Q^{Eq-^), 

and substitute q = —q^^l'^ in the coefficients of the ffist relation in ([T8|) 



-gi/2Ai(g-^E) + Ai(E) - q-^''' kx^qE) = 0, ^ Ai(^) = q-^l'^kx^qE) + q^''' Kx^q-^E). 
The second relations can be checked similarly. 

10 



The case 



The case with the harmonic oscillator is slightly complicated as the asymptotic formula 
must be modified. We utilize known facts on the Weber's function D^(z), 

(fDJz) , 1 z^ ^ 

which has an asymptotic behavior for 1] ^ 0, integer, 

Dr,{z) ^ z'^expi-zyA). 

It has the 2nd order irregular singularity at oo and regular elsewhere. 

The FSS consists of {Dr^{z), D^rj^i{iz)} or {D^(— z), D_^_i(— z^;)}. The connection 
rule reads, 



27r 

D^{z) = t^D^n-iiiz) + r''D_,_i(-2z). (19) 



There exist recurrence relations, 

D'^{z) = z/2D^{z) - D,+iiz) = -z/2Dr,{z) + vD^-i{z). (20) 

Obviously, y{x, E) is given in terms of Dfj{z). In order to cancel the phase factor arising 
from the asymptotic behavior, we define precisely 

„, ^k/2+k/2Ek ^Vk ( V 23^9 ) 

where Ek = Eq"^^, 2rjk + 1 = Ek and tj = tjq. They constitute our FSS. 

By definition, T{ {E)yQ = yi + y-i- Remembering q = i so that rji = rj_i = —t] — 1 , 
we rewrite this into the form, 



-tz) 



^^D,{z) = i^D_,_,{iz) + r^D_i_,( 
with z = y/2x. By comparing this with (|l^) we conclude 

^ ^ ' r(r/ + l) r(E/2 + l/2) 

which coincides with the result from CFT0, Q. The expectation values of Q±{E) are 
proportional to Weber's function and its derivative at the origin. Thanks to the recursion 
relations (pO]) , we can replace the latter by again Weber's function with the unit shift in 

Q_(E)ocD^(0), Q+(E)ocD^+i(0). 
11 



We shall utilize the following integral representation for Dr){z), 

D(z) = _ r(^ + ^) e--V4 /'g-tV2-.t(_i)-(.+i)^i^ 
2vrz Jc 

where C surrounds the positive real axis counterclockwise. The evaluation at z = is then 
straightforward, 

'^ > r((l-r^)/2) T{{3-E)/2y 
Hence, 



For general values of i, the representations of T| or (5± by elementary functions 
are not known. Still, we can evaluate them, e.g., from solutions to TEA ([T3|) . To fix one 
missing parameter r there, we next consider the spectral problem. 

5 Spectral Determinants and Stokes multipliers 

The final section for the n = 1 case is devoted to the spectral problem for i = 2M and 
M being an integer. We will still use E instead of A. 

We first put some remarks on elementary facts. Let H{x) be a our Hamiltonian oper- 
ator, i7 (a;) = -^ + x^^. 

Definition 1. We call ^(x) the eigen-function and E, the eigenvalue of H{x) if 
H{x)i/je{x) = EipE^x) and iPe{x) is a vector in the Hilbert space satisfying, e.g., 
\\ip{x)\\ < oo. 

Definition 2. Let P be a spatial inversion operator such that Pf{x) = f{—x) for any 
operators or vectors. 

Obviously, [H{x), P] = 0. Thus ii H{x)iPe{x) = EtpE^x) then PipE^x) = pipE^x). Since 
P^ is an identity operator, PipE^x) = iPe{—x) = ±iPe{x). Consequentially, we have 

Lemma 1. If H{x)i)E{x) = EiPe{x) then ipEix = 0) = or ^^^^U=o = 0. 

The above lemma does not require the boundary condition lim^i^ioo |V'-E(a^)| = im- 
posed by our potential. 

Two conditions can not be satisfied simultaneously. Or otherwise, ipE (^ = 0) = for 
arbitrary n, resulting a trivial tp. Thus a lemma follows. 

Lemma 2. Eigenvalues are classified by the parities of the associated eigenf unctions. We 

denote Ej' if — ^ — |a;=o = and Ej if ipE'i^ = 0) = 

12 



On the positive real axis, ipEix) = yo + a{E)yi up to normalizations. We have 
[ipE{x),yo] = a{E) from the obvious asymptotic behavior, \im.x^oo4'E{x)/yi = a{E). 
As the eigenfunction must be bounded, a{E) = ii E E {-^/} U {E~}. Conversely, if 
a{E') = for some E' then ipE'ix) is proportional to yo. Thus ipE'{x) is bounded as 
X -^ +00 and it is recessive as x ^ — cxd due to the parity argument. In addition, it 
is a solution to the eigenvalue equation. Then, by definition, E' belongs to the set of 
eigenvalues. We conclude. 

Lemma 3. // ipEi^) = Vo + ci{E)yi on the positive real axis, then a{E) oc D{E) : = 
D4E)D4E) where D±{E) := [1^(1 - E/Ef). 

Finally we quote results from the WKB analysis. 

Lemma 4. For the potential x^*^, the energy levels Ek and the spectral determinant D{E) 
behave asymptotically as 

ho{EkY ~ 27r(A; + l/2), k ^ 00 (22) 

\nD{E) ^ 77-^,E^ (23) 

2 sm(/i7r) 



M + 1 , 7ri/2r( 



1 



We shall apply the above general observation to results obtained in the preceding 
sections. The connection rule enables a representation of D{E) in terms of yo- To check 
this, we consider the Stokes matrix Sq .It connects FSS on the positive and the 
negative real axes. We start from the negative real axis. If E takes an eigenvalue, then 
'4'e{x) = yu+i apart from a normalization. The connection rule demands it behave on the 
positive axis, 

ll){x, E^) ri}^_i(g^^a)l/0 + 1"m (^a)l/l- 

Lemma ^ tells tj^,/ oc D(E). 

On the other hand, we consider the j = M case of (0). Note that g*-'^+^ = —1 , and 
yM+i = wi—x,E). Then t)^/ = i{y{x, E)y'{—x, E) — y{—x, E)y'{x, E)). Since the Ihs is 
independent of x and the rhs is not singular at a; = 0, we conveniently put x = in the 
rhs and find r|^ oc y{0, E)y'{0, E)\x=o. Thus, as a function of E, |/ol/oU=o has only zeros at 
eigenvalues. Then the above lemma leads to their identification with D^{E). The choice 
of the evaluation at x = here and in the previous section is now clear. 

Summarizing, we have a theorem. 

Theorem 5. The fusion hierarchy contains D{E) as its M— th member, 
riJ^ oc D{E), equivalently TJi\E) oc D{-E). 

A base of FSS and its derivative at the origin are proportional to spectral determinants 
depending on parities, 

y{0,E)^D4E), y'{0,E)<xD+{E). 
13 
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Figure 1: the contourplot for M = 3, \e^/'^T^^\E)\ 

The previous exphcit resuh (^) for M = 1(£ = 2) is quite consistent with this. 
Q±{E)^s are nothing but D±{E) here. They are vanishing at known spectra of the 
harmonic oscillator, E = 2n + 1, where n = (even/odd) corresponds to the parity = 
(even/odd). 

These identifications lead to the expression for T]; via spectral determinants, 



Ti'\E) = q 



1/2 



D+{Eq') , ./,D+{Eq 



D+{E) 



+ q 



-1/2 ^-W 



D+iE) 



D-{E) 



,1/2 



D_{Eq-^ 
~^JE) 



More significantly, we have BAEs, 



De{E,q^ 
D,{E^q- 



-r 



±. 



These equations, combined with the WKB result, are efficient enough to determine the 
spectral determinants, being transformed into coupled nonlinear integral equations. We, 



however, take a different route here and utilize them as a tool to investigate TEA (0). 

Let us revisit to the assumption |^ raised in section ^. Suppose all energy levels are 
enumerated exactly so that D±{E) are constructed. Then Tf '{E) is estimated. By the 
use of the analogue of the relation(|^), we can successively generate TJ (E) and check 
the validity of the assumption. Strictly speaking, as we have infinitely many levels, this 
procedure can not be accomplished. One however knows that the WKB approximation 
is fairy accurate for higher energy levels. Thus we input first 100 exact energy levels and 
approximate rests by the WKB results, to evaluate D±{E). 

Our numerical results indicate the remarkable patterns. 

Conjecture 1. Zeros of T^ (E) are of the first order and always distribute on the neg- 
ative real E axis. 

This supports the assumption, although by no means a proof. As an example, the 
contourplot for M = 3, je^/^Tg {E)\ is depicted in Fig.[l|. These patterns imply that the 
state corresponds to "the vacuum" (the ground state) in IM [|^, Q. 
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There is also an independent support to this conjecture 0, ^, ^. As shown in 0|, 
zeros of T- (E) coincide with negative of eigenvalues associated to "PT-symmetric Hamil- 
tonian p^ + x'^^{ixy with e = 2M — 2j. The numerical and analytical studies on the PT- 



symmetric Hamiltonian in [|31], ^ conclude positive and real eigenvalues for M > 1, which 
is consistent with the conjecture. The studies reveal, at the same time, the breakdown of 
the conjecture, when M being continued to a real number less than 1 O, RTl R2 . 



We assume that the validity of integral equations (|T3|) . Then all we have to do is 
fix r in evaluating Stokes multipliers and spectral determinants. The theorem ^ tells 



TJi^\-Ek) = 0, which resuhs YM^-Ekq^) = -1 or log YmI^+^I ) = {2k + l)7ii. Remember 
E = exp{9/fi). For large values of 9, numerical data concludes that the contribution form 
the integral is negligible so that we have an approximation, log yM(^) ~ '"^m exp(6'). Thus 
T)J{—Ej) = means mMexp(6'fc) = tumEJ^ = {2k + 1)tc for large enough j. Comparing 
this with the WKB result, we conclude m^^ = &o, which derives the desired quantity. 
Summarizing the results for 2nd order ODE, we have the following correspondence, 



energy <^=^ spectral parameter 
Stokes multipliers, D{E) <^=^ fusion transfer matrices 

y\x=o,y'\x=o -^^ (vacuum) expectation values of Q± 



In the next two sections, the higher order ODE will be briefly discussed. Our results 
indicate a natural extension of the ODE/IM correspondence examined for n = 1 above. 

6 Functional relations in Stokes multipliers 

As in the case of n = 1, we can introduce generalized Stokes matrices connecting dis- 
joint sectors for arbitrary n. The obvious recursion relation leads to functional relations, 
however, among complex objects corresponding to Young tableaux of the hook shape. 
Then the restriction of relations among Young tableaux of the rectangular shape results 
the desired relation ^j. This procedure requires some technique in integrable models e.g., 
quantum analogue of the Jacobi-Trudi formula. 

We derive the same relation in a simpler way using the Wronskian representation of 
Stokes multipliers. Let auxiliary functions r™ (A) be 



'TmW = [yi,y2,---ya-l,y0,ya+m,ya+m+l---yn+m\- (25) 

Note that we adopt the abbreviation defined in section 2. 

Due to yn+i+e = (— )"'?/o, Tm (A) = for m > i+ 1, £ + 2, ■ ■ ■ . This is an analogue to the 
quantum group reduction. Remark that the set contains the original Stokes multipliers as 
m = 1 cases. 

Then the claim is the following functional relations , 
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Theorem 6. 



-(0) 



(«^) _ (' 1 \a-l 



(26) 



where t^"' = lyT^""' = (—1, 

We utilize a lemma in |Q for the proof of the above theorem. 
Lemma 5. 

[fi, /2, ■ ■ ■ In, ao, fli] [/i, /2, ■ ■ ■ In, «2, 0-3] — 
[/i, /2, ■ ■ ■ /Af, ^0, 02] [/i, /2, ■ ■ ■ /at, ai, 0.3] + 
[/i, /2, ■ ■ ■ /tv, ao, %] [/i, /2, ■ ■ ■ /tv, «i, ^2] = 0. 



(27) 
(28) 
(29) 



When A^ = 1, this relation follows from the Laplace expansion of the trivial relation, 



= det 



0.2 a3\ 





/' Oq a!^ a'2 ttg 

/" a"o a"i a" 2 a" 3 

ai / 02 as 

a'l /' a'2 a^j 

\0 



a"i /" a" 2 a" 3/ 
With the similar argument, the validity of the above lemma is verified for arbitrary A^. 

Proof, of theorem ^ : 

We shall adopt identifications. 



(/i,/2, ■ ■ ■ , In) 
(ao, ai, a2, a^) 



(1/1, ■■ ■ , Va-l, Va+m+l, 
(l/O) Va, Va+m,) 2/n+m+l)- 



) Vn+r, 



For a = 1, the left hand side of the first relation should read as (1/25 
elements in eq.(^) are interpreted as 



-(«+i). 



Tt\X), (-l)"-Vr(Ag) 

(-l)"-'^riti(A), (-l)'^-Vi"l,(Ag) 



\n-2^(a-l)^ 



l)-V^"-^^(Ag), 



y^+i). The six 



D 



respectively. This immediately leads to the theorem. 

We have a note. The substitutions r™ (A) -^ (— l)"^^Tm (Ag^""*""^^/^) brings the equa- 
tion to the form known as the T— system for An in solvable models. 



Tt\M'^')TJ;:\\q-'/') = Tt^'\X)Ti^-'\X)+T!^UX)T}:U>^). 



(30) 
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This observation supports the ODE/IM correspondence for n arbitrary. In view of the 
solvable models, Tm (A) should be understood as the (eigenvalues of) transfer matrix 
associated to the auxiliary space Wm {X)- ( As the module in classical Lie algebra, Wm {X) 
is isomorphic to mA^, of which Young diagram takes a rectangular shape.) 

The relation also finds a connection to a discrete soliton system. We parameterize 
A = g^/^ and denote /(p, a, m) = Tm (X). Let Di, i = 1, 2, 3 be Hirota operators acting on 
i th variable. Then the eq (^) reads 

(exp Di - exp D2 - exp D^)/ ■ / = 0. 



This equation is known as the Hirota-Miwa equation with Zi = —Z2 = —Z^ = 1 
The present construction imposes the periodicity and boundary conditions. 



/(£ + n + 1, a, m) 
/(p, -l,m) 



/(0,a,m) 

f{p,n + 2,m) = f{p,n + l,i+l) 
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7 DVF for arbitrary n 

The DVF in the Stokes multipliers are also found for arbitrary n. 

Let us check this for r} (A) . The following lemma is useful for this purpose. 

Lemma 6. For m >2 we have a recursion relation among ratios of determinants, 

[yo, ^2, • • • Vm] ^ [yo,y2,---ym-i] [yo,yi,- ■ ■ym-i][y2,- ■ ■ ,ym] 
[yi, ■■■ ,ym] [1/1, • • • , ym-i] [yi,--- ,ym][yi,--- , ym-i] 

We should interpret [yi] -^ yi and [yo, yi,- ■ ■ Vm-i] -^ Vo for m = 2. 
Proof. The lemma is equivalent to 



(31) 



bo,i/2,---i/, 



m [yi, 



z/m-i] = [yo,y2,---ym-i][yi,--- ,ym] + [yo,yi,---ym~i] 



ii 1 ymj 



(32) 



Since ( P^ is linear in |/q , it suffices to show the equality of the coefficients of them 
in the both sides. First consider the coefficient of yo, We need to show the equality, 

[yi,--- ,ym-i][y2,---y'm] = [y2,--- ,ym][y[,---y'^^i] + [yi,--- ,ym][y2,- ■ -yln-i]- (33) 
To verify this, we prepare a matrix Ai 



M 



( z/i, y2 
y'l, 1/2 



(m— 1) (m— 1) 

\yi 1/2 



ym \ 
ym 



ym / 
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rows and ji,j2, 



Denote by 

[jl, J2, ■■ 

the minor, the determinant of a matrix obtained by deleting ii,i2 
columns from A^. Then eq. (|33|) is represented as, 



D 



where D = detA^. Obviously this is the Jacobi identity. Thus the equality of coefficients 
of yo in both sides is established. The equalities are similarly proven up to those of i/q . 



m 
m 


D 


'1 
1 


= D 


'1' 
m 


D 


m 
1 


+ DD 


"1, 
.1, 


m 
m 



For t/q case, the first term of the rhs in (|32D does not contribute. We can check 



however, the equality of the reminding terms. Thus the lemma is proved. 

Next we will show 
Theorem 7. r{ ^(A) can be represented in the following DVF 

-(i)mn _ [y2,y3,--- ,yn+i] , [y2,y3,--- ,yn][yo,yi,--- ,yn-i\ 



D 



rr^(A) 



[yi, 1/2, ■ ■ ■ , yn] [yi, 1/2, ■ ■ ■ , yn-i] [yi, z/2, 
[y2,y3][yo,yi,y2] , [yo,yi]y2 , yo 



Z/nJ 



+ 



+ 



[1/1,1/2] [1/1,1/2,1/3] [1/1, 1/2] 1/1 1/1' 



(34) 



Proof. Firstly we substitute r{ 

/l/l, 1/2, 

y'l 1/2, 



(n+l) 



^— l)" to eq.(pD and obtain 



yn \ 

y'n 



\yt^ yt\ 



,{n) 



rf'(A) 



yrj \rr{x)/ 



(n), 



/i/o-(-iri/„+i\ 

i/^-(-i)X+i 

W^ - (-i)"2/S J 



(1) 



The application of Cramer's formula yields t{ in the form 



-(1) 



[l/2,l/3,--- ,1/n+l] [1/0,1/2, 



.Vn 



[l/l,l/2,--- ,1/n] [1/1,1/2, ■■■ ^VnY 

We use Lemma ^ to the second term in the rhs to obtain, 

(1) _ [l/2,l/3,--- ,1/n+l] , [1/0, l/l,--- ,l/n-l][l/2,l/3,--- ,1/n-l] , [l/O, 1/2, - - - , 1/n-l] 



Ti 



[1/1,1/2,--- ,1/n] [1/1,1/2, --- ,1/n] [1/1,1/2,--- ,1/n-l] [l/l, 1/2, 



1/n-lJ 



It is now obvious that repeated applications of Lemma ^ to the last term results the 
expression (PH). D 
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We mimic the case of n = 1 and introduce D functions 
Then the a + 1— th term in (|^) reads, 

[y2,l/3,---,l/a+l]b0,l/l,---,l/a] _ a-n/2 D^''^'\\q''/')D^''\\q^''^'^/') 

The DVF consists of n + 1 terms for the solvable Uq{A\ ) model of which auxiliary 
space is Ai as a classical module. It is characterized by Baxter's Q operators of n species, 

^(1).,^ ^ g(")(Ag("+^)/^) g(n)(Ag(n-3)/2)g(.-l)(^g./2) 

Q('^+i)(Ag('^-2)/2)QW(Ag('^+i)/2) 

+J- g(a+l)(Ag'^/2)g(a)(Ag{a-l)/2) + ' ' ' 

Q(2)(^^-i/2)g(i)(^^) g(^)(Ag~^) 

^^' g(2)(Agi/2)g(i)(Ag) + ^° g{i)(A) ' 

Clearly, we have 

Theorem 8. Under the identification, r{ (A) ^^ T{ (Ag), 

gW(A)^D(")(A), /a^g-"/=^+" (35) 

two DVFs coincide. 

The pole-free property of r} (A), required from the linear independence of FSS, results 
BAE, 



Q(«-i)(A5")gi/2)g(a)(A(.-)g-i)g(a+i)(AWgi/2)' ' ' ' ^ 

where Q^^^iXf^) = and Q("+i) = g(°) = 1. 

Thus we have verified the common algebraic structure for arbitrary n. 

The representation (0) or the identification ( pSj) is , however, not unique. One easily 
recognizes this by remembering the simplest case {n = 1) where two different expressions 
are available for rj . This originates from the simple fact that both y and its derivative 
are solutions to Baxter's T — Q relation. The situation is also true for n > 1. One can 
show that the identification 

gW(Ag-+^-) ^ £-[^,^,,^,^2, . . . ,|/fc+J|,=o, J = 0, 1, • • • (36) 
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works and we have a variety of representations for the same r} . This is shown by using 
formulas analogous to Lemma ^ and the detail will be published elsewhere. 

Before closing the section, we comment on Tm , {a > 1). The corresponding DVFs are 
known in the integrable models, but explicit forms are quite involved. Still, one can param- 
eterize them by analogue of Young tableaux||28||. We utilize the "tableaux" representation 



in proving the equivalence of the DVF in integrable models and Stokes multipliers Tm for 
a, n and m general. This point may be further discussed in a separate publication. 

8 Summary and discussions 

In the present report, we discuss a curious connection between n + 1 the order ODE 
and integrable models. When n = 1, the connection is efficient enough to derive analytic 
equations which yields estimations of eigenvalues and Stokes multipliers. 

For higher n, the correspondence is still at the algebraic relation level. Unfortunately, 
the definition of the eigenvalue problem is not necessary clear for higher order differential 
equations. The characterization of the eigenspace (it is the Hilbert space for n = 1) is not 
obvious. More technically, there are several subdominant solutions in each sector. This 
obscures the identification of D{E) in the general Stokes matrices. The lack of the con- 
nection prevents us from writing down the integral equations and evaluating parameters 
like "r" for n = 1. We however comment some progress made in [^ ^ 

The observation made in the last few sections may be interesting. Suppose that the 
ODE/IM correspondence even occurs at the construction of models. Then one may find 
the variable x also in IM. Once if one of Baxter's Q^^^ is constructed, the other indepen- 
dent Q^^' functions are found in derivatives of the Q^^' with respect to the hidden variable 
X. Moreover, one can generate higher Q functions, i.e., Q^°'\a > 2 mere by taking deter- 
minants of fundamental Q^^\ The other Q'^'^^'s are again obtainable via taking derivatives. 
On the other hand, construction of Q functions via the standard "pair-propagation" ar- 
gument fl^ seems to be far more complicated for ra > 1. To the authors' knowledge, the 



explicit construction of Q is done only for cases corresponding n = \ [|T^, |3^, ^, and the 
procedure is already involved. The systematic construction found in ODE is not obvious 
in IM. The present results for ODE may be a clear guide for analyses in the analogous 
issue in IM, but it needs further research. 

Finally, we comment that the ODE/IM correspondence is still at the " phenomenolog- 
ical" stage. The fundamental question as to the origin of the correspondence is still open. 
The complete classification of the ODE tractable with the IM approach may need the 
answer to this fundamental question. 

I hope to clarify these issues in future publications. 
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